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Abstract 

What determines the caliber of axonal branches? We pursue the 
hypothesis that the axonal caliber has evolved to minimize signal 
propagation delays, while keeping arbor volume to a minimum. We 
show that for a general cost function the optimal diameters of 

mother (c/q) and daughter {^d^,d^ branches at a bifurcation obey 

a branching law: d^*^ = d^*^ + Jj^^ . The derivation relies on the 
fact that the conduction speed scales with the axon diameter to the 
power V ( V = 1 for myelinated axons and V = 0.5 for non- 
myelinated axons). We test the branching law on the available 
experimental data and find a reasonable agreement. 

1 Introduction 

Multi-cellular organisms have solved the problem of efficient transport of nutrients 
and communication between their body parts by evolving spectacular networks: 
trees, blood vessels, bronchs, and neuronal arbors. These networks consist of 
segments bifurcating into thinner and thinner branches. Understanding of branching 
in transport networks has been advanced through the application of the optimization 
theory ([1], [2] and references therein). Here we apply the optimization theory to 
explain the caliber of branching segments in communication networks, i.e. neuronal 
axons. 

Axons in different organisms vary in caliber from 0.1 |im (terminal segments in 
neocortex) to lOOOjxm (squid giant axon) [3]. What factors could be responsible for 
such variation in axon caliber? According to the experimental data [4] and cable 
theory [5], thicker axons conduct action potential faster, leading to shorter reaction 
times and, perhaps, quicker thinking. This increases evolutionary fitness or, 
equivalently, reduces costs associated with conduction delays. So, why not make all 
the axons infinitely thick? It is likely that thick axons are evolutionary costly 
because they require large amount of cytoplasm and occupy valuable space [6], [7]. 
Then, is there an optimal axon caliber, which minimizes the combined cost of 
conduction delays and volume? 



In this paper we derive an expression for the optimal axon diameter, which 
minimizes the combined cost of conduction delay and volume. Although the relative 
cost of delay and volume is unknown, we use this expression to derive a law 
describing segment caliber of branching axons with no free parameters. We test this 
law on the published anatomical data and find a satisfactory agreement. 

2 Derivation of the branching law 

Although our theory holds for a rather general class of cost functions (see Methods), 
we start, for the sake of simplicity, by deriving the branching law in a special case 
of a linear cost function. Detrimental contribution to fitness, C, of an axonal 
segment of length, L , can be represented as the sum of two terms, one proportional 
to the conduction delay along the segment, T , and the other - to the segment 
volume, V : 

C = aT + /JV. (1) 

Here, (X and are unknown but constant coefficients which reflect the relative 

contribution to the fitness cost of the signal propagation delay and the axonal 
volume. 
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Figure 1 : Fitness cost of a myelinated axonal segment as a function of its diameter. 
The lines show the volume cost, the delay cost, and the total cost. Notice that the 
total cost has a minimum. Diameter and cost values are normalized to their 
respective optimal values. 

We look for the axon caliber d that minimizes the cost function . To do this, we 
rewrite C as a function of d by noticing the following relations: i) Volume, 
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V= — Ld ; ii) Time delay, T= — ; Hi) Conduction velocity S = kd for 

4 s 
myelinated axons (for non-myelinated axons, see Methods): 



C = a- + j3L-d=L 
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This cost function contains two terms, which have opposite dependence on d , and 
has a minimum. Fig. 1 . 

Next, by setting = U we find that the cost is minimized by the following axonal 

dd 

caliber: 
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(3) 



The utility of this result may seem rather limited because the relative cost of time 
delays vs. volume, (xj P ,is unknown. 




Figure 2: A simple axonal arbor with a single branch point and three axonal 
segments. Segment diameters are do, dj, and d2. Time delays along each segment are 
to, t], and t2. The total time delay down the first branch is Ti= to +t], and the second - 

T2= to +t2. 



However, we can apply this result to axonal branching and arrive at a testable 

prediction about the relationship among branch diameters without knowing the 

relative cost. To do this we write the cost function for a bifurcation consisting of 
three segments. Fig. 2: 

C = ai(ta+0 + a2(ta+t2) + j3(Va+Vi+V2), (4) 

where is a conduction delay along segment 0, - conduction delay along 

segment 1, ?2' conduction delay along segment 2. Coefficients OC^ and 



represent relative costs of conduction delays for synapses located on the two 
daughter branches and may be different. We group the terms corresponding to the 
same segment together: 



(5) 



We look for segment diameters, which minimize this cost function. To do this we 
make the dependence on the diameters explicit and differentiate in respect to them. 
Because each term in Eq. (5) depends on the diameter of only one segment the 
variables separate and we arrive at expressions analogous to Eq.(3): 
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It is easy to see that these diameters satisfy the following branching law: 

(3 _ j3 , j3 



dl = d[ + d^ . 



(6) 



(7) 



Similar expression can be derived for non-myelinated axons (see Methods). In this 
case, the conduction velocity scales with the square root of segment diameter, 
resulting in a branching exponent of 2.5 . 

We note that expressions analogous to Eq. (7) have been derived for blood vessels, 
tree branching and bronchs by balancing metabolic cost of pumping viscous fluid 
and volume cost [8], [9]. Application of viscous flow to dendrites has been 
discussed in [10]. However, it is hard to see how dendrites could be conduits to 
viscous fluid if their ends are sealed. 

Rail [11] has derived a similar law for branching dendrites by postulating 
impedance matching: 



■■df^df. 



(8) 



However, the main purpose of Rail's law was to simplify calculations of dendritic 
conduction rather than to explain the actual branch caliber measurements. 

3 Comparison with experiment 

We test our branching law, Eq.(7), by comparing it with the data obtained from 
myelinated motor fibers of the cat [12], Fig. 3. Data points represent 63 branch 
points for which all three axonal calibers were available. Eq.(7) predicts that the 
data should fall on the line described by: 



+ 



d^ 



= 1, 



(9) 



where exponent T] = 3 . Despite the large spread in the data it is consistent with our 
predictions. In fact, the best fit exponent, rj = 2.57 , is closer to our prediction than 
to Rail's law, 77 = 1.5. 

We also show the histogram of the exponents TJ obtained for each of 63 branch 
points from the same data set. Fig. 4. The average exponent, Tj = 2.67 , is much 



closer to our predicted value for myelinated axons, 77 = 3, than to Rail's law, 
7 = 1.5. 
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Figure 3: Comparison of the experimental data (asterisks) [12] with theoretical 
predictions. Each axonal bifurcation (with ^ (^2) is represented in the plot twice. 

The lines correspond to Eq.(9) with various values of the exponent: the Rail's law, 
?7 = 1 .5 , the best- fit exponent, TJ = 2.57 , and our prediction for myelinated axons, 
7 = 3. 

Analysis of the experimental data reveals a large spread in the values of the 
exponent, 77. This spread may arise from the biological variability in the axon 

diameters, other factors influencing axon diameters, or measurement errors due to 
the finite resolution of light microscopy. Although we cannot distinguish between 
these causes, we performed a simulation showing that a reasonable measurement 
error is sufficient to account for the spread. 

First, based on the experimental data [12], we generate a set of diameters Jq, 

and branch points, which satisfy Eq. (7). We do this by taking all diameter 

pairs at branch point from the experimental data and calculating the value of the 
third diameter according to Eq. (7). Next we simulate the experimental data by 
adding Gaussian noise to all branch diameters, and calculate the probability 
distribution for the exponent T] resulting from this procedure. The line in Fig. 4 

shows that the spread in the histogram of branching exponent could be explained by 
Gaussian measurement error with standard deviation of 0.4/vm . This value of 
standard deviation is consistent with 0.5 /Jm precision with which diameter 
measurements are reported in [12]. 
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Figure 4: Experimentally observed spread in the branching exponent may arise from 
the measurement errors. The histogram shows the distribution of the exponent 77, 

Eq. (9), calculated for each axonal bifurcation [12]. The average exponent is 
J] = 2.67 . The line shows the simulated distribution of the exponent obtained in the 
presence of measurement errors. 

4 Conclusion 

Starting with the hypotheses that axonal arbors had been optimized in the course of 
evolution for fast signal conduction while keeping arbor volume to a minimum we 
derived a branching law that relates segment diameters at a branch point. The 
derivation was done for the cost function of a general form, and relies only on the 
known scaling of signal propagation velocity with the axonal caliber. This law is 
consistent with the available experimental data on myelinated axons. The observed 
spread in the branching exponent may be accounted for by the measurement error. 
More experimental testing is clearly desirable. 

We note that similar considerations could be applied to dendrites. There, similar to 
non-myelinated axons, time delay or attenuation of passively propagating signals 
scales as one over the square root of diameter. This leads to a branching law with 

exponent of 5/2. However, the presence of reflections from branch points and 

active conductances is likely to complicate the picture. 

5 Methods 

The detrimental contribution of an axonal arbor to the evolutionary fitness can be 
quantified by the cost, C . We postulate that the cost function, C , is a 
monotonically increasing function of the total axonal volume per neuron, V , and all 
signal propagation delays, , from soma to j -th synapse, where _/= 1,2,3, ••• : 



(10) 



Below we show that this rather general cost function (along with biophysical 
properties of axons) is minimized when axonal caliber satisfies the following 
branching law: 

rfo^=<+J2' (11) 

with branching exponent = 3 for myelinated and // = 2.5 for non-myelinated 
axons. 

Although we derive Eq. (11) for a single branch point, our theory can be trivially 
extended to more complex arbor topologies. We rewrite the cost function, C , in 
terms of volume contributions, Vj , of i -th axonal segment to the total volume of 

the axonal arbor, V , and signal propagation delay, t- , occurred along i -th axonal 
segment. The cost function reduces to: 

civ.+v.+v^jo+hA+ti)- (12) 

Next, we express volume and signal propagation delay of each segment as a 
function of segment diameter. The volume of each cylindrical segment is given by: 

v;.=^Vf, (13) 

where and rf, are segment length and diameter, correspondingly. Signal 
propagation delay, , is given by the ratio of segment length, L- , and signal speed, 
Sj . Signal speed along axonal segment, in turn, depends on its diameter as: 

Si=kd^, (14) 

where V = l for myelinated [4] and V = 0.5 for non-myelinated fibers [5]. As a 
result propagation delay along segment i is: 

t:=—^. (15) 

Substituting Eqs. (13), (15) into the cost function, Eq. (12), we find the dependence 
of the cost function on segment diameters. 
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(16) 



To find the diameters of all segments, which minimize the cost function £ , we 
calculate its partial derivatives with respect to all segment diameters and set them to 
zero: 



— = C'y-L,d,-C',^ = (17) 
dd^ 2 ' ' Mr 

By solving these equations we find the optimal segment diameters: 

, 2v(g:;+£;) , iv<t^ , 2vc; 

" yt;r£; kn€y kTrt'y 

These equations imply that the cost function is minimized when the segment 
diameters at a branch point satisfy the following expression (independent of the 
particular form of the cost function, which enters Eq. (18) through the partial 

derivatives , , and ): 

d^=d1+dl, r] = v + l. (19) 



References 

[I] Brown, J. H., West, G. B., and Santa Fe Institute (Santa Fe N.M.). (2000) Scaling in 

biology. Oxford; New York: Oxford University Press. 

[2] Weibel, E. R. (2000) Symmorphosis : on form and function in shaping life. Cambridge, 
Mass.; London: Harvard University Press. 

[3] Purves, D. (1997) Neuroscience. Sunderland, Mass.: Sinauer Associates. 

[4] Rusliton, W. A. H. (1951) A theory of the effects of fibre size in meduUated nerve. J 
Physiol 115, 101-122. 

[5] Hodgkin, A. L. (1954) A note on conduction velocity. J Physiol 125, 221-224. 

[6] Cajal, S. R. y. (1999) Texture of the Nervous System of Man and the Vertebrates, Volume 
1. New York: Springer. 

[7] Chklovskii, D. B., Schikorski, T., and Stevens, C. F. (2002) Wiring optimization in 
cortical circuits. Neuron 34, 341-347. 

[8] Murray, C. D. (1926) The physiological principle of minimum work. I. The vascular 
system and the cost of blood volume. PNAS 12, 207-214. 

[9] Murray, C. D. (1927) A relationship between circumference and weight in trees and its 
bearing on branching angles. J Gen Physiol 10, 725-729. 

[10] Cherniak, C, Changizi, M., and Kang D.W. (1999) Large-scale optimization of neuron 
arbors. Phys Rev E 59, 6001-6009. 

[II] Rail, W. (1959) Branching dendritic trees and motoneuron membrane resistivity. Exp 
Neurol 1, 491-527. 

[12] Adal, M. N., and Barker, D. (1965) Intramuscular branching of fusimotor fibers. J 
Physiol 177, 2SS-299. 



